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Abstract 

Evolutionary branching points are a paradigmatic feature of adaptive dynamics, be- 
cause they are potential starting points for adaptive diversification. The antithesis to 



in 

o 

(3) evolutionary branching points are Continuously stable strategies (CSS's), which are 

'^ convergent stable and evolutionarily stable equilibrium points of the adaptive dynam- 

rN ics and hence are thought to represent endpoints of adaptive processes. However, this 

assessment is based on situations in which the invasion fitness function determining 
the adaptive dynamics have non-zero second derivatives at a CSS. Here we show 
that the scope of evolutionary branching can increase if the invasion fitness func- 
tion vanishes to higher than first order at a CSS. Using a class of classical models 
for frequency-dependent competition, we show that if the invasion fitness vanishes 
to higher orders, a CSS may be the starting point for evolutionary branching, with 
the only additional requirement that mutant types need to reach a certain threshold 

1 



frequency, which can happen e.g. due to demographic stochasticity. Thus, when in- 
vasion fitness functions vanish to higher than first order at equilibrium points of the 
adaptive dynamics, evolutionary diversification can occur even after convergence to 
an evolutionarily stable strategy. 



Introduction 



Since the publication of the seminal papers Metz et al. ( 1992[ ), Metz et al. ( 1996| ), 
Dieckmann & Law| ( |1996[ ) and |Geritz et al.\^99%\ , the theory of adaptive dynamics 
has attracted a lot of attention as a useful tool to study long-term evolutionary dynam- 
ics in continuous phenotype spaces. In particular, the phenomenon of evolutionary 
branching emerged as paradigmatic feature of adaptive dynamics that serves as a 



basic metaphor for adaptive diversification and speciation (Dieckmann & Doebeli 



T999I IDieckmann et a/.j [20041 IDoebelij [20T0l [Doebeli & Dieckmamil [2003| [Gentz 



& Kisdi 2000 Geritz etal. 1998). During evolutionary branching, a population first 
evolves under directional selection to an evolutionary branching point in phenotype 
space, at which selection turns disruptive due to frequency dependence, which in turn 
causes the population to split into diverging phenotypic clusters. 

Evolutionary branching is a generic feature of adaptive dynamics and can be ob- 



served in a multitude of models incorporating frequency-dependent selection (Doe 



beli 2010). The analytical understanding of evolutionary branching rests on the dis- 
tinction between two basic stability concepts of adaptive dynamics, which are easiest 
to understand in the context of adaptive dynamics in 1 -dimensional trait spaces. In 
general, adaptive dynamics is derived from the invasion fitness function function 
f{x, y) describing the long-term per capita growth rate of a rare mutant type y in a 



resident population that is monomorphic for trait value x (Geritz et al. 1998 Metz 



g?a/.[|1996l). The adaptive dynamics of the trait x is then given as a gradient dynamics 



of the invasion fitness function and described by the canonical equation (Dieckmann 



&Law 1996): 



^^-M(x)^^^^'^) 



(1) 

y=x 



dt dy 

The partial derivative with respect to the mutant trait is evaluated at the resident trait, 

and the selection gradient df /dy is in general a vector of the same dimension as the 
phenotype space. M{x) is a matrix describing the mutational process in the pheno- 
type under consideration, but for 1 -dimensional trait spaces, M{x) is a scalar that 
only influences the speed of evolution. For 1 -dimensional traits one therefore often 
assumes M{x) = 1 to simplify the analysis. Equilibrium points of the adaptive dy- 
namics are then given as points x* in phenotype space at which the selection gradient 
vanishes, i.e., points x* satisfying 
df{x*,y) 



= 0. (2) 

y=x* 



dy 

A singular point is called convergence stable if it is locally stable for the adaptive 
dynamics (1), i.e., if the system (1) returns to the singular point after any perturbation 
that is small enough. On the other hand, a singular point is called evolutionarily stable 
if the invasion fitness function f{x*, y) has, as a function of y, a maximum at x* , for 
in that case no mutant trait in the vicinity of the singular point has a higher growth 
rate than the singular trait value itself. 

Even though both types of stability, convergence stability and evolutionary sta- 
bility, are derived from the invasion fitness function f(x,y), the two conditions for 
stability are in general different. For example, it is possible that a singular point 
is evolutionarily stable, but not convergence stable, a situation that is referred to as 
Garden-of-Eden configuration ( Nowak & Sigmund] |1990| ), because even though a 



resident with the singular trait value cannot be invaded by any nearby mutant, the 
resident population will actually evolve away from the singular point if it is started 
with a resident value that is arbitrarily close (but not equal) to the singular point. 



On the other hand, it is also possible that a singular point is convergent stable, 
but not evolutionarily stable, in which case the population will evolve to the singular 
point, only to find itself being vulnerable to invasion by any nearby mutant. In such 
a situation the singular point is a minimum for the invasion fitness function, and con- 
vergence to such a point therefore opens the scope for evolutionary diversification. 

When the invasion fitness function only vanishes to first order at a singular point 
(i.e., only the selection gradient, but not the higher order derivatives of the invasion 
fitness vanish at the singular point), the adaptive dynamics around a singular point 
can be fully classified in terms of the second derivative of the invasion fitness at 
the singular point ( [Diekmann] |2004[ |Doebeh[ |2010[ [Geritz et alj |1998| ). In short, a 



convergence stable singular point (and only those are really of interest) is either an 
endpoint of the evolutionary dynamics, or it is an evolutionary branching point at 
which an evolving population diversifies into diverging phenotypic clusters. In the 
first case, the convergent stable singular point is also evolutionarily stable, so that af- 
ter convergence to the singular point, the population is immune to invasion by nearby 
mutants. Such points are classically called Continuously Stable Strategies (CSS;|Ei 



hel (1983)) and represent evolutionary endpoints at least under the assumption of 
small mutations. 

In the second case, the convergent stable singular point is evolutionarily unstable, 
and it can be shown that in 1 -dimensional trait spaces, these two conditions imply 
two more conditions: that of mutual invasibility, and that of evolutionary divergence. 
Mutual invasibility refers to the fact that trait values on either side of the singular can 
mutually invade each other. More precisely, if e > is small enough, then a resident 
population that is monomorphic for trait value x* — t can be invaded by a mutant 
x*+e, and vice versa. In particular, phenotypes on either side of the singular point can 
coexist in a protected polymorphism. Moreover, in such a polymorphic population 
consisting of traits x* — e and x* + e, evolution in each monomorphic subpopulation 



occurs away from the singular point, so that the phenotypic distance between the two 
coexisting clusters increases, and hence the clusters diverge evolutionarily. 

When the invasion fitness function only vanishes to first order at a singular point, 
a CSS may satisfy the condition for mutual invasibility, but it never satisfies the con- 
dition for evolutionary divergence of coexisting subpopulations. Thus, while pheno- 
types on either side of a CSS may coexist, evolution in two coexisting phenotypic 
clusters always induces convergence of both clusters to the singular trait value. In 
particular, a CSS cannot be the starting point for evolutionary diversification. 

These statements can be made precise in terms of the second derivatives of the 
invasion fitness function at the singular point ( [Diekmann , 2004 Geritz et al. 1998 1. 
In particular, the condition for convergence stability of the singular point, i.e., for 
local stability of an equilibrium of the dynamical system (1), is given by 



d 

dx 



df{x,y) 



dy 



y=x=x' 



d^f{x,y) 



dxdy 



d'^f{x,y) 



y=x=x' 



dy^ 



<0. 



(3) 



y=x=x' 



(4) 



Also, the condition for evolutionary stability of a singular point is simply 

d^f{x,y) ^Q 

dy'^ 

a y=x=X* 

If both these conditions are satisfied, the singular point x* is a CSS, and hence no di- 
versification occurs. On the other hand, if the first equality is satisfied but the second 
is reversed (so that the singular point is a minimum of the invasion fitness function), 
then not only is the singular point convergent stable, but both mutual invasibility and 
evolutionary divergence as described above are satisfied, and hence x* is an evolu- 
tionary branching point (Diekmann, 2004[ Doebeli, 2010| ). 

In this paper, we show that the scope of evolutionary branching can increase if the 
invasion fitness function vanishes to higher than first order at a singular point. More 
precisely, using a class of classical models for frequency-dependent competition we 
show that if the invasion fitness vanishes to higher orders, a CSS may satisfy both 
the conditions for mutual invasibility and the condition for evolutionary divergence. 



and hence can be the starting point for evolutionary branching. This can happen 
even if mutations are small, with the only additional requirement that mutant types 
need to reach a certain threshold frequency, e.g. due to demographic stochasticity. 
Thus, when invasion fitness functions vanish to higher than first order at equilibrium 
points of the adaptive dynamics, evolutionary diversification can occur even after 
convergence to an evolutionarily stable singular point. 



Model and Results 

We consider models for frequency-dependent competition in which individuals 
are characterized by a 1-dimensional phenotype x (e.g., body size). The phenotype 
determines two ecological properties. On the one hand, the carrying capacity of 
populations that are monomorphic for trait value x is given by the carrying capacity 
function K{x), and on the other hand the strength of competition between two in- 
dividuals with phenotypes x and y is given by the competition kernel a{x,y). Note 
that by definition, the carrying capacity is a property of populations, but the fact that 
K{x) is the carrying capacity of monomorphic populations makes it a function of 
the phenotype x (moreover, it is easy to define this function at the individual level in 
terms of per capita birth and death rates; see e.g. Chapter 3 in Doebeli (2010)). 



The ecological dynamics of a resident population that is monomorphic for phe- 
notype X is assumed to be logistic: 

riVj 1 - — ^ , (5) 



dt "" \ K{x 

where A^a; is the population size of the resident x, and r > is the intrinsic growth 

rate, which is assumed to be independent of the phenotype x. Without loss of gener- 
ality we set r = 1 in the following. Clearly, A^^: = K{x) at ecological equilibrium. 
When a mutant y appears in a resident x, the dynamics of its population density A^j^ 



is given by 



''^'-jv,fi^^' + °!'';">^-|. (6) 



* ' V K(y) 

where the competition kernel a{y^ x) is used to describe the competitive impact of 
the resident x on the mutant y. If the resident is at its ecological equilibrium K{x) 
and the mutant is rare, so that A^j^ is negligible, the per capita growth rate of the 
mutant becomes 

/(a;, y) is the invasion fitness function, which determines the adaptive dynamics. 

Here we consider the following class of functions for the carrying capacity and 
the competition kernel: 

A-W =AVxp (-1 (M)') (8) 

{ \ { \x — y\\'^\ 
a(x,y)=exp(^--^— — J J (9) 

We assume that the exponent is a real number s > 2 to ensure that K is at least twice 
differentiable at and a is at least twice differentiable aty = x. (Note that the first 
derivative of fi' at and of a at y = x is 0.) The carrying K(x) decreases with 
increasing distance from its unique maximum at a; = at a rate that is determined 
by the parameter a^- In evolutionary terms, K{x) describes a stabilizing component 
of selection whose strength is measured by ax (with small ax indicating strong 
stabilizing selection). The parameter Kq scales the population density. 

The competition kernel a{x, y) describes the frequency-dependent component of 
selection, and for any given trait value x, the strength of the competitive impact of 
individuals with trait value y decreases with increasing distance \y — x\ at a rate that 
is determined by the parameter a a, which therefore determines the strength of fre- 
quency dependence (with small a^ indicating a rapid decline of competitive impacts 
with phenotypic distance). Note that the competition kernel is symmetric in the sense 



that it only depends on the phenotypic distance \y — x\. Also note that for s = 2, 
the carrying capacity and the competition kernel have Gaussian form, which results 
in a classical model for frequency-dependent competition that has a long tradition 
in ecological and evolutionary theory (Dieckmann & Doebeli 1999[ [KTrkpatrick & 
Barton! |7997| [MacArthur & Levinsj [T9671 IRoughgardenj [T9791 ). 

With the invasion fitness function (7), the adaptive dynamics (1) becomes 
dx df{x,y) 

y=x 



dt 



dy 
da{y,x 



dy 



+ 



y=x 



K\x) 
K{x) 



(10) 



K'{x) 
' K{x) 



(where for simplicity and without loss of generality we have assumed M{x) = 1 
for the mutational process). Expression (10) reflects the fact that away from singular 
points, the adaptive dynamics is solely determined by the stabilizing component of 
selection, i.e., by the carrying capacity function K{x). This is essentially due to the 
assumption that the competition kernel is symmetric and differentiable at y = x, 
which implies that its derivative at x = y must be 0. It immediately follows that 
x* = is the only singular point of the adaptive dynamics, because x* = is the 
only solution of K'{x*) = 0. Moreover, even though the second derivative of K{x) 
at X* = may be 0, and hence the Jacobian of the dynamical system (10) at the 
equilibrium x* = may be 0, it is clear that x* = is a globally stable equilibrium 
for the dynamics (10), because dx/dt > for x < x* and dx/dt < for x > x*. 
Thus, X* = is a convergent stable singular point. 

Evolutionary stability of the singular point is determined by the second derivative 
of the invasion fitness function. For s = 2, we have 



d'f{0,y) 



dy^ 



crt. 



a 



(11) 



y=0 ^a ^K 

Therefore, x* = is evolutionarily stable (i.e., the partial derivative (11) is negative). 



and hence a CSS, if and only if a^ > ctr, which is a familiar result for the Gaussian 



case s = 2 ( Dieckmann & Doebeli[ [T999| ). Also, x* is an evolutionary branching 



point, and hence the starting point of adaptive diversification, if and only if a a < ctk 
The situation changes, however, for s > 2. 
For s > 2, we have 



= 0, (12) 

y=o 



but it is easy to see that, again, the invasion fitness f{x*, y) has a maximum aty = 
if and only if a a > (Tk- Thus, the evolutionary stability condition is the same 
for all s > 2. In particular, the singular point is a CSS if and only if a^ > (Tk- 
What changes, however, is that for s > 2, the conditions for mutual invasibility and 
for evolutionary divergence that are necessary for evolutionary branching may be 
satisfied even if aa > ctr, i-C-, even if the singular point is a CSS. 

To see this, we assume cTq, > ax, and we first consider coexistence of two strate- 
gies e > and — e close to the singular point x* = 0. (We note that we make 
the symmetry assumption of considering phenotypes e and — e for analytical conve- 
nience, and that the arguments given below are also valid for resident strains e and 
5, as long as the two strains lie on either side of the singular point and |e| and |5| are 
small.) The per capita growth rate of a rare mutant — e in a resident e at equilibrium 
K{e) is 

a(—e,e)K(e) 
1 - V / = 1 - "(-^' ^) > 0' (13) 

and similarly for the growth rate of a rare e mutant in a — e resident. It follows that 
the two types can coexist, and if N^ denotes the equilibrium population density of 
the two subpopulations at the coexistence equilibrium (the two subpopulations have 
the same density at equilibrium due to symmetry), one easily calculates that 



K {e){l-a{e,-e)) _ K{e) 
1 — a(e, — e)^ 1 + a(e, — e) 



N. = 7^ _r^\,^^ = . : r . ■ (14) 



To check for evolutionary divergence, we assume two coexisting residents e and — e 
and calculate the invasion fitness /(e, — e, y) for rare mutants y occurring in either 



resident: 



/(e, -e,l/) 



N,{a{y,e) + a{y,-e)) 
K{y) 



(15) 



A salient example of this invasion fitness function is shown in Figure 1, which illus- 
trates that selection in the two resident strains e and — e points away from the singular 
point, i.e., that selection is divergent. 

This can be made mathematically precise by considering the selection gradients 
in the two resident strains. For example, it is fairly easy to show that 



dfie,-e,y) 



dy 



se 



y=e 



-^\s] 



In particular, this selection gradient vanishes for e* satisfying 

2 (l + exp[2'~\^y] 



(16) 



a 



(17) 



K 



an equation that can be solved to yield 



.s-l 



aK 



(Tr 



e* = 2 — a^ i In 
Clearly, a solution e* > exists if 



3-2 

(Tq < 2 = aK- 



-1 



(18) 



(19) 



It follows that if inequality (19) is satisfied, there will be some resident pair (e, — e) 
with e small enough for which the selection gradient in the resident e will be posi- 
tive, and hence selection will drive this resident away from the singular point (since 
we assume e > 0). By symmetry, the selection gradient in the resident — e will be 
negative, so that overall, selection will drive evolutionary divergence of the two coex- 
isting strains. Moreover, (e*, — e*) is the equilibrium of the adaptive dynamics of the 
two coexisting and diverging strains, where e* is given by (18). At this equilibrium, 
both resident strains are maxima for the invasion fitness function (15), as illustrated 
in Figure 1 . 
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Figure 1a 



Figure 1b 



Figure 1c 




mutant trait z 



mutant trait z 



E 1.1 

mutant trait z 



Figure 1 : Invasion fitness as a function of the mutant trait z when the resident is monomorphic 
for the singular strategy x* (panel a), and when the resident population consists of two coexisting 
strategies e and — e (panels b,c). Panel a shows that the singular point is a fitness maximum (and 
hence a CSS, see main text). Panel b shows that in the invasion fitness function with two residents e 
and — e, the maximum at x* =0 shown in panel a resolves into two maxima, one lying to the right 
of e and the other lying to the left of — e. As a consequence, the derivative of the invasion fitness 
is positive at e and negative at — e, leading to divergent evolution (note that the fitness function is 
necessarily at the two resident strategies). Panel c illustrates that divergent evolution comes to halt 
when the residents come to lie on fitness maxima themselves, i.e., when they reach the values e* and 
— e* given by eq. (18). Parameter values: s = 4, ct/^ = 1, (Ta = 1.2, Kq = 1. 
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Overall, assuming s > 2 in the ecological functions (8) and (9), we now have the 
following situation. For ax < (Ja < 2'^^^'^'>/^aK, the singular point x* is a CSS, i.e., is 
convergent stable and evolutionarily stable, but two strains e and — e with e > small 
enough can mutually invade each other and hence coexist, and selection generates 
evolutionary divergence in such coexisting strains. Thus, for a^ < o-^ < 2^^^'^'>/^aK 
the singular point is both a CSS and satisfies the conditions for evolutionary branch- 
ing. Figure 2 summarizes these results. 



Figure 2 



CSS without branching 




Evolutionary branching 



4 
exponent s 



Figure 2: Parameter regions for evolutionary branching for adaptive dynamics based on the eco- 
logical functions (8) and (9). For s > 2 there is a region of parameter space, bordered by the curves 
Ca/cK ~ 1 and OajoK — s^"^^'/", in which the adaptive dynamics (10) has a CSS at x* =0 that 
can be the starting point of evolutionary branching. The point indicated by a cross corresponds to the 
parameter values used for Figures 3 and 4. 

However, there is a snag, which is that while phenotypes on either side of the sin- 
gular may invade each other for gk < cr^, they cannot invade the singular strategy, 
which is, after all, evolutionarily stable under these conditions. As a consequence, the 
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establishment of a coexisting coalition (e, — e) cannot occur deterministically from 
very low frequencies. However, establishment of such coalitions can occur when 
the density of the mutants e and — e reaches a certain threshold. In fact, the ecolog- 
ical systems consisting of three strains e, (the singular strategy) and — e exhibits 
bistability, as can be seen as follows. 

Given three strains (e,0, — e), and assuming an idealized invasion scenario in 
which both e and — e have density 5 and the singular strain has density K{0) — 26 = 
Kq — 26, the per capita growth rate of the marginal strains is given by 

^ a{e, 0){Ko - 26) + 6{a{e, e) + a{e, -e)) 



K{e) 
while the per capita growth rate of the singular strain is 



(20) 



^ _ a(0, 0)(i^o - 2^) + 6{a{0, e) + ajO, -e)) ^^^^ 

Because the singular point is an ESS (always assuming ax < cto), it follows that 
(21) is larger than (20) for 6 (and e) small enough, so that the ecological equilibrium 
at which only the singular strain has positive density is locally stable. However, as 
6 becomes larger and reaches a certain threshold, the growth rate of the marginal 
strains becomes larger than the growth rate of the singular strains. By equating the 
two growth rates (20) and (21), one can find an analytical expression for the threshold 
value of 5*, 

S* = Ko^^^^p^^ + 0{e^) (22) 

and it is possible to show that for densities of the marginal strain above this thresh- 
old, the growth rate of the marginal strains is always larger than the growth rate of 
the singular strain. This implies that once the density of the marginal strains has 
crossed the threshold, the singular strain will go extinct, and coexistence of the coali- 
tion (e, — e) will be established. Thus, the ecological equilibrium at which only the 
marginal strains have positive density is also locally stable. 
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In fact, this threshold value can be fairly small, in which case one would expect 
that it can be reached due to random fluctuations in models incorporating demo- 
graphic stochasticity. Indeed, it is easy to observe evolutionary branching after con- 
vergence to a CSS in individual-based implementations of the competition models 
described here. In such models, individuals are characterized by their trait value, and 
their death rate is determined by the carrying capacity on the one hand, and by the 
effective density on the other hand. For an individual with trait value x, the effective 
density Nejj is obtained by summing over all other individuals in the population, 
where each individual y is counted with a weight a(x, y): N^ff = J2y '^(^i v)- The 
death rate of individual x is then Neff/K{x), and all individual are assumed to have 
a birth rate of 1 . Individual birth and death rates then define a stochastic evolutionary 
process, as described in detail in Dieckmann et al. ( 2004| ) (Chapter 4) and in Doe 



beli ( 2010[ ). When an individual gives birth during this process, the phenotype of the 



offspring is chosen from the mutation kernel, a normal distribution that has a mean 
equal to the parent value and some small variance amut- Note that in such a model, 
all individuals present at any point in time have different phenotypes, and when the 
average phenotype of the population is close to the singular value, the population 
actually consists of a cluster of phenotypes around the singular value. Because birth 
and death events are stochastic, phenotypes at a distance > e from the singular point 
may reach non-negligible densities due to random fluctuations even if they are not 
favored deterministically by selection. As a consequence of these stochastic fluctu- 
ations, evolutionary branching can occur after convergence to a CSS in accordance 
with the analytical arguments given above. This is illustrated in Figure 3. 
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Figure 3 
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Figure 3: Evolutionary branching from a CSS in the individual-based model. Parameter values 
were s — 4, ctk = 1, (Tq = 1.2, Kq ~ 270, and <Tmut = 0.01 for the width of the mutation kernel. 
The population was initialized with individuals drawn from a Gaussian distribution with mean —1 and 
a small variance. 
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Discussion 



The theory of adaptive dynamics (Dieckmann & Law 1996 Geritz et al. 1998 



Metz et al.\ 1996) has established itself as an excellent mathematical toolbox for 



understanding long-term evolutionary dynamics of quantitative traits. In particu- 
lar, the paradigmatic phenomenon of evolutionary branching has been useful for un- 
derstanding the ecological conditions leading to adaptive diversification and specia- 
tion dPieckmann & Doebelil \T999\ [Dieckmann "gForj [20041 poebelij [20T0l |Kisdi & 



Geritz 1999). During the process of evolutionary branching, a population evolves 
under directional selection to a point in phenotype space where selection turns dis- 
ruptive due to frequency-dependence, after which the population splits into diverging 
phenotypic clusters. 

The point at which selection ceases to be directional is a so-called singular point, 
and if the trait under consideration is 1 -dimensional, i.e., given by a real number, 
the conditions for evolutionary branching can be described by properties of the inva- 
sion fitness function at the singular point. The first derivative of the invasion fitness 
function vanishes by definition at a singular point, and if the second derivative is 
non-zero, then the singular point is an evolutionary branching point if and only if 
it is convergent stable and evolutionarily unstable ( Diekmami] 2004, DoebeU[ 2010[ 



Geritz etaL\ [1998[ ). These conditions already reveal that evolutionary branching is 
a generic phenomenon that occurs not only at particular points in parameter space, 
but for whole parameter regions. Here we showed that the realm of evolutionary 
branching can be even increased if the invasion fitness vanishes to higher than first 
order in the Taylor expansion around a singular point. More precisely, under these 
conditions a Continuously Stable Strategy (CSS), i.e., a convergent stable and evo- 
lutionarily stable singular point, can be the starting point for evolutionary branching, 
which is never the case when the second derivative of the invasion fitness function at 
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the singular point is non-zero. 

The mechanism for the observed phenomenon can be appreciated graphically 
by considering a perturbation of the invasion fitness function at the singular point 
(Figure 1), i.e., by comparing the two invasion fitness functions F{z) = f{x* — 
e, X* + e, z) and F{z) = f{x*, z), where x* is the singular point, e > is small, and 
z is the phenotype of a rare mutant. If x* is a CSS, then F has a maximum at x* . If 
dF{x*)/dz 7^ 0, then to second order, F is a parabola with maximum at x*, and it 
follows that the perturbed function F(z) is also a parabola with negative curvature, 
and with zeroes at x* — e and x*+e. This fitness function therefore generates selection 
towards the singular point in both perturbed resident strains. 

However, if dF{x*)/dz = 0, then the multiplicity of the maximum at x* is higher 
than 1, and as a consequence, the perturbed function F(z) has two maxima in the 
vicinity of x*, one lying to the left of x* — e and the other lying to the right of x* + e. 
Therefore, selection in the two perturbed residents is away from the singular point, 
as illustrated in Figure lb. Thus, in terms of selection, when the second derivative of 
F at X* is zero resolution of the multiplicity of the maximum of F at x* can result in 
a scenario that is opposite to the situation occurring when the second derivative of F 
at X* is non-zero (and hence x* is a "regular" CSS). 

We have described the parameter range for which this happens using a general 
class of models for frequency-dependent competition. In particular, CSS's, which 
are usually thought to be endpoints of the evolutionary dynamics, can be the start- 
ing point for evolutionary diversification due to frequency-dependent competition. 
Recently Frans Jacobs has made similar observations for invasion fitness functions 
that are given as fourth-order polynomials around the singular point ( [Jacobs ????). 



In our models, the exponent of the first non-zero order in the Taylor expansion of 
the invasion fitness function around the singular is the same as the exponent s ap- 
pearing in the carrying capacity function (8) and in the competition kernel (9). We 
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have assumed that s > 2 to ensure twice differentiability of the invasion fitness at 
the singular point. For s = 2, the model converts to the familiar Gaussian model, in 
which the singular point is a CSS for ax < cf^, and an evolutionary branching point 
for ax > cTa, where ax and a a are the curvatures of the competition kernel and the 
carrying capacity. Because the second derivative of the invasion fitness function at 
the singular point is non-zero for s = 2, no evolutionary diversification occurs in this 
case when the singular point is a CSS, i.e., when ax < o^- For s > 2, the singular 
point is still a CSS for gk < o"a, but because the second derivative at the singular 
point is 0, evolutionary branching can occur for ax < Ua < 2'--^^'^^/^aK, i.e., for a 
range of parameters for which the CSS condition holds. Note that for all s evolution- 
ary branching still occurs for ax > cTa- (In fact, one can show that for exponents 
< s < 2, the threshold ax = era remains the boundary between evolutionary di- 
versification on the one hand, and CSS's as endpoints of the evolutionary dynamics 
on the other hand.) Thus, compared to models in which the invasion fitness func- 
tion only vanishes to first order at the singular point, the region of parameter space 
leading to evolutionary diversification is enlarged when the invasion fitness function 
also vanishes to second order. We note that this phenomenon also occurs when the 
exponents in the carrying capacity and the competition kernel are different. Based 
on both analytical and numerical studies, and denoting these exponents by sx and 
Sa, respectively, we conjecture that branching can occur after converging to a CSS 
whenever Sq, > 2 and sx > Sa- 

To conclude, we mention an interesting connection of the adaptive dynamics 
models studied here to the corresponding logistic partial differential equations mod- 
els 

-^r = ^^(") [' im — ) ' ^'^^ 

which describe the dynamics of phenotype distributions (f){x) due to frequency-dependent 
competition. Here K(x) and a{x,y) are the carrying capacity function (8) and the 
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competition kernel (9), and the integral J a{x, y)(l){y)dy is the effective density expe- 
rienced by individuals of phenotype x. As before, the dynamics (23) can be studied 
for different exponents s of the ecological functions (8) and (9). For the Gaussian 
case s = 2, it is well known that system (23) has an equilibrium distribution con- 
sisting of a single delta peak at x* = for a^ < o"q, and a Gaussian equilibrium 
distribution with maximum at a;* = and a positive variance for ax < era- This 
corresponds to the results from adaptive dynamics in the sense that the boundary 
for aK = Ua separates two different dynamic regimes, one in which diversity is 
not maintained and another in which it is, even though maintenance of diversity 
in this latter regime occurs in the form of unimodal phenotype distributions, rather 
than in the form of multimodal distributions representing distinct phenotypic clus- 
ters. [Pigolotti et aL\ ( |2009| ) have shown that the Gaussian case s = 2 is structurally 
unstable, and that for s > 2 the regime in which diversity is maintained is charac- 
terized by multimodal equilibrium distributions. Here we can add that not only does 
model (23) have equilibrium distributions consisting of multiple delta peaks when 
s > 2, but such diversification can occur even for cta' < era, as illustrated in Figure 4. 
In fact, our numerical simulations indicate that diversification can occur in the partial 
differential equation model (23) exactly for the same range of parameters as in the 
adaptive dynamics model, i.e., for aa < 2^^^'^^/^aK- 

Moreover, for ax < <Ja < 2^^^'^'>/^aK, and in complete agreement with the adap- 
tive dynamics model, the dynamics (23) converges to a single delta peak at x = if 
the initial phenotype distribution is very narrow, and diversification can only occur if 
the initial distribution is sufficiently wide. Overall, however, the partial differential 
equation models confirm that for exponents s > 2 in the ecological functions, the 
region of parameter space leading to adaptive diversification is enlarged. 



19 



Figure 4 



5000 




ecological trait x 

Figure 4: Evolutionary diversification in the form of a bimodal equifibrium distribution branching 
in the logistic partial differential equation model (23). Parameter values were the same as for Figure 
3. A convolution with a Gaussian mutation kernel of width 0.01 was incorporated in the birth term in 
(23) to make the temporal evolution consistent with the individual-based simulation shown in Figure 
3. Note that the position of the two modes at equilibrium corresponds well with the two phenotypic 
clusters resulting from evolutionary branching in Figure 3. 
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